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IMAGE PROCESSING METHOD, APPARATUS, AND STORAGE MEDIUM 
FOR RECOGNITION OF IRRADIATION AREA 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

The present invention relates to an image 

? processing method, ai/ image processing apparatus, and a 
storage medium using secondary differences for 
recognition of an/ irradiation field. 
10 Related Background Art 

The recent progress in the digital technology 
makes it common to convert a radiographic image to a 
digital image signal, subject this digital image signal 
to image processing, and display it on a CRT or the 
15 like or print it out. In photography of the 

radiographic image it is common practice to diaphragm 
the irradiation area so as to radiate the radiation 
only into a necessary area of a human body on the 
humane grounds or for the reason of preventing 
20 scattering from unnecessary areas so as to prevent 

lowering of contrast. Normally, a processing parameter 1 
is determined from distribution of density values of 
the image prior to the image processing and the image 
processing is then carried out based on the parameter 
25 thus determined. However, if the irradiation area is 
not limited, unwanted information, in a sense, outside 
an area of interest is used for the determination of 
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the image processing parameter, which poses a problem 
that appropriate image processing is not effected. 

It is thus necessary to extract the irradiation 
area and determine the image processing parameter from 
5 information of only the area of, interest. As an 

extraction method of the irradiation area, for example, 
Japanese Patent Application (Laid-Open) No. 05-007579 
describes a method for dividing the image area into 
small areas and extracting the irradiation area, based 

10 on values of variance in the small areas. Further, for 
example, Japanese Patent Publication No. 6-90412 
describes a method for expressing a change of image 
densities between a predetermined number of pixels near 
an end portion of the image area by an approximate 

15 equation substantially of a linear equation and 
extracting the irradiation area, based on the 
difference between assumed image density values 
according to this approximate equation and actual 
density values indicated by a sample image signal. 

20 These methods are predicated on the premise that 

an object area in the photographic image includes the 
irradiation area (an area actually irradiated), and a 
process of judging whether the object area is an area 
including the irradiation area (which will also be 

25 referred to as an "area with an irradiation diaphragm" ) 
or an area not including the irradiation area (which 
will also referred to as an "area without the 
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irradiation diaphragm" ) is carried out as a preliminary 
operation before execution of these methods. 

Methods for judging presence/absence of the 
irradiation diaphragm include methods for comparing an 
5 average, a median, or the like of densities in a 

central area of the image with an average of densities 
of the object area and judging that the object area is 
an area with the irradiation diaphragm, if the average 
of densities of the object area is not more than a 
10 predetermined value, as described in U.S. Patent No. 
5,091,970, for example. 

SUMMARY OF THE INVENTION 

The method of above Japanese Patent Application 

15 (Laid-Open) No. 05-007579 had to obtain the density 

variance in each small area and thus had the problem of 
increased computation complexity and computation time. 
Even within the area of interest, density values vary 
suddenly, for example, at the edge of the lang field 

20 and a change rate of densities could become higher 
there than at the edge of the irradiation area. 
Particularly, the variances become high in the area 
where the lang is in contact with the ribs at the edge 
of the lang. This caused candidate points for the 

25 irradiation area to be extracted from outside the edge 
of the irradiation area as well and this posed another 
problem that the judgment became difficult. 



In the method of .above Japanese Patent Publication 
No. 6-90412, a foot area outside the irradiation area 
was assumed for calculation of the linear approximate 
equation and there arose the problem that the method 
was invalid if this foot area was not able to be 
extracted well. Further, it was also assumed that the 
foot area outside the irradiation area included two 
areas of a gentle portion and a quickly rising portion 
and this posed another problem that the method was 
invalid where the entire foot area was gentle or where 
there were three or more changes of density gradient. 

The present invention has been accomplished in 
order to solve the problems described above and an 
object of the invention is, therefore, to extract the 
irradiation area accurately. 

An image processing method of the present 
invention is an image processing method comprising: 

a step of determining a plurality of areas 
arranged in a predetermined direction on an image and 
each having a predetermined shape; 

a step of calculating a secondary difference value 
of density values representing the respective areas in 
the plurality of areas; and 

a step of judging one end point of an irradiation 
area from the secondary difference values calculated in 
the calculating step. 

In the conventional image judgment method as 
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described in above U.S. Patent No. 5,091,970 etc., for 
example, where the photographic image is one in a state 
in which the subject overlaps with the edge of the 
image, the average density value at the end of the 
5 image varies depending upon the area of the subject 

over the edge of the image and upon the transmittance 
of the radiation. Therefore, the presence/absence of 
the irradiation diaphragm was sometimes misjudged in 
the area at the end of the image. 

10 When the radiant dose of the irradiation diaphragm 

was small, the density difference also became small 
between the central part of the image and the object 
area and the presence/absence of the irradiation 
diaphragm in the object area was misjudged in some 

15 cases. 

The present invention has been accomplished to 
eliminate the drawbacks described above and an object 
of the invention is, therefore, to make it possible to 
accurately judge whether an image area as a processed 

20 object is an area including the irradiation area. 

An image processing method of the present 
invention is an image processing method for extracting 
an irradiation area in an input image, the image 
processing method comprising: 

25 a step of detecting an irradiation end, based on a 

density distribution in each area, for a plurality of 
areas in a desired direction in the image; and 
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a step of evaluating the result of the detection, 
based on the result of irradiation ends detected for 
each of the plurality of areas. 

Another image processing method of the present 
invention is an image processing method for judging 
whether an object area in an image includes an 
irradiation area, the method comprising: 

a secondary difference vaJ?ue acquisition step of 
acquiring secondary difference values from one- 
dimensional image data of the object area; 

an irradiation end extraction step of extracting a 
coordinate of an end point of/ the irradiation area from 
the secondary difference va/ues acquired in the 
secondary difference valuer acquisition step; 

a comparison step of comparing the coordinate 
extracted in the irradiation end extraction step with a 
coordinate of an end point of the irradiation area 
included in the image, the coordinate being obtained 
preliminarily; and 

a judgment step of judging whether the object area 
includes the irradiation area, based on the result of 
the comparison in the comparison step. 

The other objects, arrangements, and effects of 
the present invention will become more apparent with 
the description of the embodiments of the present 
invention which follows. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Fig, 1 is a block diagram of an irradiation area 
extraction device according to Embodiment 1; 

Fig. 2 is a block diagram to show the structure 
5 inside an irradiation area extraction unit; 

Fig. 3 is a flow chart to show processing 
procedures in the irradiation area extraction unit; 

Fig. 4 is a structural diagram to show calculation 
areas; 

10 Fig. 5A arid Fig. 5B are structural diagrams to 

show a radiographic image and density values of one 
line on the radiographic image; 

Fig. 6A and Fig. 6B are characteristic diagrams to 
show density values of one line on Fig. 5A and values 
15 of secondary difference on this line; 

Fig. 7 is a block diagram to show the structure of 
an image judgment device of Embodiment 2-1; 

Fig. 8 is a diagram for explaining an example of 
image data inputted into the above image judgment 
20 device; 

Fig. 9 is a flow chart for explaining a processing 
program carried out in the above image judgment device; 

Fig. 10 is a block diagram to show the structure 
of an image judgment device of Embodiment 2-2; 
25 Fig. 11 is a flow chart for explaining a 

processing program carried out in the above image 
judgment device; 
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Fig* 12 is a block diagram to show "the structure 
of an image judgment device of Embodiment 3-1; 

Fig. 13 is a diagram for explaining an example of 
image data inputted into the above image judgment 
5 device; 

Fig. 14 is a flow chart for explaining a 
processing program carried out in the above image 
judgment device; 

Fig. 15 is a block diagram to show the structure 
10 of an angle extraction device and an area extraction 
device in Embodiment 4-1; 

Fig. 16 is a flow chart of* a processing procedure 
sequence in the angle extraction device and the area 
extraction device; 
15 Fig. 17 is a diagram for explaining the processing 

in the angle extraction device and the area extraction 
device; 

Fig. 18 is a flow chart of a processing procedure 
sequence for obtaining an approximate center in an area 
20 in Embodiment 4-4; 

Fig. 19 is a block diagram to show the structure 
of an area extraction device in Embodiment 4-7; 

Fig. 20 is a flow chart of a processing procedure 
sequence in the area extraction device of Embodiment 4- 
25 7; 

Fig. 21 is a block diagram to show the structure 
of an area extraction device in Embodiment 4-9; and 
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Fig. 22 is a flow chart of a processing procedure 
sequence in the area extraction device of Embodiment 4- 
9. 



5 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

( Embodiment 1 ) 

Fig. 1 is a block diagram to show the structure of 
an irradiation area extraction device according to an 
embodiment of the present invention. In the figure, 

10 reference numeral 101 designates an image input unit 
for accepting input of an image from the outside, 102 
an irradiation area extraction unit for extracting an 
irradiation area irradiated with radiation from the 
image inputted at the image input unit 101, 103 an 

15 image processing unit for performing image processing 
based on the irradiation area extracted at the 
irradiation area extraction unit 102, and 104 a display 
unit for displaying an image subjected to the 
processing at the image processing unit 103. 

20 Fig. 2 is a block diagram to show the inside 

structure of the irradiation area extraction unit 102, 
in which reference numeral 200 designates a calculation 
area input unit for accepting input of a direction, a 
start point, and an end point yfor determination of 

25 calculation areas and 201 a calculation area 

determination unit for determining, based on the input 
information at the calculation area input unit 200, 
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areas of calculation carried out by a/calculation unit 
202 from the input image supplied fr?bm the image input 
unit 101. Numeral 202 denotes a calculation unit for 
calculating primary difference values and a secondary 
difference value, described hereinafter, from the 
calculation areas determined aft the calculation area 
determination unit 201, 204 A judgment unit for judging 
irradiation area ends, 203 /a memory unit for storing 
the values calculated at yche calculation unit 202 and 
10 the irradiation area encas judged at the judgment unit 

204, and 205 an irradiation area determination unit for 
determining the irradiation area from the irradiation 
area ends stored in the memory unit 203 and judged at 
the judgment unit: 204. 
15 Fig. 3 is a flow chart to show fhe flow of the 

processing at the irradiation area ^extraction unit 102, 
Fig. 4 is a diagram to show the calculation areas 
determined at the calculation ar^ea determination unit 
201. Fig. 5A is a diagram to afhow a radiographic image 
20 and Fig. 5B is a diagram to show density values on line 
X 0 to X 3 of Fig. 5A. The abscissa indicates positions 
on the line X 0 to X 3 and thp ordinate density values on 
the line. Fig. 6A is an Enlarged view of part of Fig, 
5B. Fig. 6B is a plot of secondary difference values 
25 calculated at the calculation unit 202 against points 
on the above line between X 0 and X 3 . 

The operation will be described below. 
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In Fig. 1, an image from an external device not 
illustrated is inputted through the image input unit 
101 into the irradiation area extraction unit 102. The 
irradiation area extraction unit 102 extracts the 
irradiation area from the input image and delivers the 
information of the irradiation area to the image 
processing unit 103, The image processing unit 103 
performs the image processing based on the information 
of the irradiation area extracted at the irradiation 
area extraction unit 102 and the image through the 
image processing is displayed on the display unit 104. 

The processing at the irradiation area extraction 
unit 102 of Fig. 2 will be described next according to 
the flow of Fig. 3. The irradiation area extraction 
unit 102 receives the input of the calculation 
direction, the calculation start point, and the 
calculation end point necessary for /the determination 
of the calculation areas at the calculation area input 
unit 200 (step S301 ) . It is, however, noted that the 
input herein is required only fcyr an initial input 
screen and it does not have to yr>e set for the next 
input screen and after. Next/ the calculation area 
determination unit 201 determines the calculation 
areas, based on the information from the calculation 
area input unit 200. The/calculation areas herein are 
three areas A, B, and C /of a rectangular shape arranged 
in parallel on the imajge, as illustrated in Fig. 4, and 
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are used as calculation areas for calculation of the 
secondary difference value at a calculation point 
indicated by (x,y). These three areas A, B, and C have 
an equal area which is deterirfined by two parameters of 
5 "a" and "b" and the distance between the areas is 

represented by "d" . These parameters "a" , "b" , "d" are 
determined experimentally . The calculation area 
determination unit 2t)l preliminarily determines the 
calculation areas/ for all calculation points and the 
10 result, together with information of duplicate 

calculation/areas, is stored in the memory unit 203 
( step S3p2 ) . 

Then the calculation unit 202 calculates 
representatives of density values im the above 
15 calculation areas A, B, and C according to a method 
described hereinafter and further /calculates the 
primary difference values and ther secondary difference 
value. The primary difference values are two values, 
which are a value "e" resulting/ from subtraction of the 
20 density representative in the area B from the density 
representative in the area C and a value 11 f" resulting 
from subtraction of the density representative in the 
area A from the density representative in the area B. 
The secondary difference value is a value resulting 
25 from subtraction of the value "f" from the value "e" . 

The primary difference values and secondary difference 
value calculated in this /way are stored in the memory 
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unit 203. From the duplicate infc&rmation of the 
calculation areas stored in the /memory unit 203 , the 
primary difference values for duplicate calculation 
areas are obtained from the values stored in the memory 
5 unit 203 and the secondary difference value thereof is 
calculated using them ( steps S303 and S304 ) . 

Then the judgment unit 204 judges thfe irradiation 
area ends from the primary difference values and the 
secondary difference values stored in/the memory unit 
10 203. Fig. 6B shows the secondary difference values 

calculated at the calculation unlit 202. The secondary 
difference values are negatively large at the points X 1 
and X 2 where the densities vary suddenly, and the 
aforementioned primary difference values f take 
15 positive values in a density increasing direction (when 
seen from the side of the/point X 0 ) but take negative 
values in a density decreasing direction. From this 
property the judgment unit 204 judges that a point 
where the above secondary difference value is minimum 
20 and the primary difference values are positive is a 

candidate for an ^hd point of the irradiation area. If 
there are plural /candidates, the judgment unit 
determines that? a first appearing candidate is an end 
point of the /irradiation area ( step S305 ) . Then the 
25 irradiation/ area end point is stored in the memory unit 
203 (step/S306). 

Further, the judgment unit ^04 judges whether the 
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processing is completed for all the directions supplied 
from the calculation area input unit fZOO and, if not, 
the processing is repeated from step S302 (step S307). 
After the irradiation area end points are determined in 
all the directions, the irradiati'on area determination 
unit 205 determines the irradiation area. Here, a line 
passing the end point of the irradiation area and being 
perpendicular to the calculation direction of the 
secondary difference value por example, the direction 
from X 0 to X 3 in Figs. 5A and 5B) is calculated for all 
the end points of the irmdiation area and an area 
surrounded by these lines obtained is determined as an 
irradiation area ( stejy S308 ) . 

The present embodiment presents the effects of 
less computational complexity and/capability of readily 
extracting the end point of the /irradiation area, • 
because the calculation unit Z02 is arranged to 
calculate the representative^ of density values in the 
calculation areas A, B, C amd calculate the secondary 
difference value based opt the representative values. 

Since the device of the present embodiment has the 
irradiation area determination unit 205 for determining 
the irradiation area, the present embodiment has such 
an effect that the irradiation area can be extracted 
even if the irradiation area is rectangular, polygonal, 
or circular. The present embodiment further has such 
an effect that when the irradiation area is polygonal, 
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an almost circumscribed area can be extracted by- 
selecting directions of lower dimensions than it. 

Described below are methods for calculating the 
representative value in each of the above calculation 
5 areas A, B, C. Simple methods for such calculation 
include a method for calculating an average of all 
density values (all pixel values) in a calculation 
area, a method for sorting all the density values in a 
calculation area according to the density values and 

10 using a density value of a pixel at an intermediate 
position, i.e., a median, and so on. There are also 
methods for using the limited number of pixels in each 
calculation area, instead of all the pixels, by 
sampling or the like according to the necessity and 

15 using an average of density values of the limited 

number of pixels or a median of the limited number of 
pixels . 

Here, when the method for calculating the 
representative of density values in the calculation 

20 area is the method using the average of all the density 
values in the calculation area, the method has the 
effect of being resistant to noise. When the 
representative is determined by the method for sorting 
all the density values in the calculation area 

25 according to the density values and using the density 
value of the pixel at the intermediate position, the 
method has the effects of simple calculation and being 
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resistant to noise. Further, when the representative 
is determined by the method using the average of the 
limited number of points in the calculation area or the 
method for sorting the density values of the limited 
5 number of points and using the density value at the 
intermediate position as a representative value, the 
methods have the effect of capability of further 
decreasing the computational complexity* 

Another approach for determining a representative 
10 of each calculation area A, B, C shown in Fig. 4 will 
be described below. First, let us assume that there 
are the calculation areas having the vertical height of 
"a" and the horizontal height of "b" , the distance 
between the areas is represented by "d", and each 
15 calculation point is indicated by (x,y), as illustrated 
in Fig. 4. 

In the present embodiment, where the 
representative value of eafch area A, B, or C is S(A), 
S(B), or S(C), respectively, the secondary difference 
20 value SS(X) is calculated according to Eq. (1) below. 

SS(X) = S(A) - 2 x S(B) + S(C) ... (1) 

Next described is a method for calculating the 
representative value according to this approach where 
the representative of each area is indicated by S(X) in 
25 general. First, a pixel value at each calculation 

point (x,y) is defined as f(x,y) and values Fl(x) are 
calculated according to Eq. (2) below. 
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Fl (x) = f*f (x,y)dy 
J o 



(2) 



Here, the integration range is the vertical height 
"a" of the above rectangular region and for simplicity 
of the description thereof Eq. (2) is expressed so as 
5 to perform the integration in the range of "0" to "a". 
Further, values F2(x) are next calculated 
according to Eq. (3) below. 

F2(x) = min{Fl(x + xl ) - h(xl)|xieK} 

... (3) 

10 Then S(x) is calculated using the values F2(x), 

according to Eq. (4) below. 

S(X) = max{F2(x - xl ) + h(xl)|xl€K} 

... (4) 

Here, h(x) is a function represented by Eq. (5) 
15 below and K is a domain of definition thereof. 
h(x) = 0 , -b/4 < x < b/4 

= -oo / otherwise ... ( 5 ) 

Then each of the left and right irradiation ends 
is obtained from a value of x taking minSS(X) in the 
20 left half and the right half of the image. The upper 
and lower ends can also be obtained in a similar 
fashion. 

In other words, the above approach is nothing but 
an operation in which projection of pixel values is 
25 made with respect to a- predetermined direction (the 

vertical direction) of /the area, values obtained by the 
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projection are smoothed using a one-dimensional 
morphology filter, and yche secondary difference value 
is calculated with -tl\d distance d. 

With this calculation method, the irradiation area 
5 can be detected with less errors than in the case of 
the above methods using the average and median. 

The functional blocks 101 to 104, 200 to 205 in 
Fig. 1 and Fig. 2 may be configured on a hardware basis 
or may be configured in a microcomputer system 
10 comprised of CPU, memory, and so on. When they are 
configured in the microcomputer system, the above 
memory comprises the memory medium according to the 
present invention and this memory medium stores a 
program for carrying out the processing illustrated in 
15 the flow chart of Fig. 3. This memory medium can be 

selected from semiconductor memories such as ROM, RAM, 
and the like, optical disks, magneto-optical disks, 
magnetic media, and so on, or it may be substantiated 
in the form of either one of a CD-ROM, a floppy disk, a 
20 magnetic tape, a non-volatile memory card, and so on. 
Since the device according to the present 
embodiment is arranged to determine the calculation 
areas comprised of the plural areas of the 
predetermined shape arranged in the predetermined 
25 direction, calculate the secondary difference value of 
the density values representing the respective areas in 
the plural areas, and judge one end point of the 
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irradiation area from the secondary difference value 
thus calculated, the present embodiment presents the 
effects of capability of decreasing the computational 
complexity and in turn ^decreasing the computation time 
5 and capability of extracting the end points of the 
irradiation area wyrh accuracy. 

Since the device of the present embodiment is 
arranged to determine the irradiation area from the 
plural end points of the irradiation area judged, it 
10 can enjoy the effect of capability of extracting the 
irradiation area even if the irradiation area is 
rectangular, polygonal, or circular. When the 
irradiation area is polygonal, the device presents the 
effect of capability of extracting the almost 
15 circumscribed area by selecting directions of lower 
dimensions than it. 

Since the density value representing each area in 
the plural areas is an average density value in each 
area, the device can enjoy the effect of capability of 
20 extracting the end points of the irradiation area with 
accuracy even from an image with noise. 
( Embodiment 2 - 1 ) 

The present embodiment is carried out, for 
example, by an image judgment device 2100 illustrated 
25 in Fig. 7. This image judgment device 2100 is provided 
with a judgment unit 2110, a control unit 2120 for 
carrying out control of the operation of the judgment 
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unit 2110, and a program memory 2130 to which the 
control unit 2120 makes access, as illustrated in Fig. 
7. 

The judgment unit 2110 is composed of /a secondary 
difference value calculation circuit 2201/ for 
calculating the secondary difference value from data of 
an object area in an input image, a left end point 
extraction circuit 2202 for extracting a left end point 
of an irradiation area included in /the object area, 
10 based on the secondary difference/value calculated at 
the secondary difference value calculation circuit 
2201, a right end point extract/ion circuit 2203 for 
extracting a right end point of the irradiation area 
included in the object area,/ based on the secondary 
15 difference value calculated^ at the secondary difference 
value calculation circuit /2201 , and a diaphragm 
presence/ absence judgment circuit 2204 for judging 
whether the object area/ is an area with an irradiation 
diaphragm or an area without an irradiation diaphragm, 
20 from the left end poi/nt extracted at the left end point 
extraction circuit 2202 and the right end point 
extracted at the r£ght end point extraction circuit 
2203. 

Here, for example, an image G as illustrated in 
25 Fig. 8 is supplied as an input image to the judgment 
unit 2110. This input image G is a two-dimensional, 
radiographic, thoracic part image obtained by 
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photographing the -thoracic part with the irradiation 
diaphragm by use of an imaging device with the 
irradiation diaphragm function. 

In above Fig. 8, "X a " and "X b " represent positions 
of end points (irradiation area ends) of the 
irradiation area preliminarily obtained with respect to 
the horizontal axis H, and "A" to "D" are areas at the 
edges of the image. 

These areas A to D are object areas for which 
whether the irradiation diaphragm is present or absent 
(i.e., whether the irradiation area is present or 
absent) is judged. In this example, the areas A and B 
are areas with the irradiation diaphragm while the 
areas C and D are areas without the irradiation 
diaphragm. 

Various processing programs for controlling the 
operation of the judgment unit 2110 are preliminarily 
stored in the program memory 2130. 

Specifically, for example, a processing program 
according to the flow chart as illustrated in Fig. 9 is 
preliminarily stored in the program memory 2130 and 
this processing program is read and executed by the 
control unit 2120, so as to operate the judgment unit 
2110 as follows. 

First, when the area C is an object area, the 
secondary difference value calculation circuit 2201 
calculates the secondary difference values SS(x) from 





... (21) 

In this Eq. (21), " f(x)" indicates image data of a 
one-dimensional line which crosses the area C in the 
horizontal direction, "x" coordinate values thereof, 
and "d" a constant indicating a difference distance. 

Next, using the secondary difference values SS(x) 
calculated at the secondary difference value 
calculation circuit 220l/ the left end point extraction 
circuit 2202 extracts jche left end point xl of the 
irradiation area included in the area C according to 
Eq. (22) below (step S2302 ) . 

SS(xl) = min{SS (x) t 0 <> x <; x} (22) 

x; coordinates on the horizontal axis between the 

left and right ends of the object area (coordinates 
between xa and xb ) 

On the other hand, using the /secondary difference 
values SS(x) calculated at the secondary difference 
value calculation circuit 2201 J the right end point 
extraction circuit 2203 extracts the right end point x2 
of the irradiation area included in the area C 
according to Eq. (23) bel/w (step S2303 ) . 

SS(x2) = min{SS(x) Ix <. x <. Length} 

... (23) 

In this Eq. (23), "Length" indicates the length of 

£3 
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the area C along the horizontal axis. 

Then the diaphragm presence/absence judgment 
circuit 2204 obtains orders of approximation of the 
left end point xl extracted at the left end point 
5 extraction circuit 2202 and the right end point x2 
extracted at the right end point extraction circuit 
2203 with respect to the left and right irradiation end 
points preliminarily obtained (the positions X a and X b 
of the irradiation ends along the horizontal axis H) 

10 and, if the orders of approximation are high, the 

diaphragm presence/absence judgment circuit 2204 judges 
that the left end point xl and the right end point x2 
are the left and right end points of the irradiation 
area (the area irradiated directly), i.e., that the 

15 area C is an area without the irradiation diaphragm. 

On the other hand, if the orders of approximation are 
low, it is judged that the area C is an area with the 
irradiation diaphragm (step S2304). 

Specifically, using a constant "e" indicating an 

20 approximation width, if the below conditions are 

satisfied the area C is judged as an area without the 
irradiation diaphragm; otherwise the area C is judged 
as an area with the irradiation diaphragm, 
-e < xl - X a < e and 

25 -e < x2 - X b < e 

When in an object area the irradiation ends are 
detected similarly to the irradiation ends in the 
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adjacent area, the object area is judged as an area 
irradiated with radiation as the adjacent area was. 

Steps S2301 to S2304 are carried out for each of 
the other areas A, B, D, in a similar fashion to the 
area C. 

As described above, the c/resent embodiment is 
arranged to calculate the secondary difference values 
SS(x) from the data of the object area subjected to the 
judgment of presence /absence of the radiation diaphragm 
in the two-dimensional input image G and judge whether 
the object area/is an area with or without the 
irradiation cJ^Laphragm, using the secondary difference 
values SS^k) . 

Construction of the device in this structure 
permits boundary points between the directly irradiated 
area and the other areas to be extracted with accuracy 
even from photographic image obtained when a subject 
with low transmittances of radiation is photographed. 
Therefore, presence or absence of the irradiation 
diaphragm can be judged with accuracy for the object 
areas including the irradiation area. The presence or 
absence of the irradiation diaphragm can also be judged 
for an object area in a photographic image in which a 
portion with low radiation transmittances, such as the 
abdominal part or the like, overlaps with the edge of 
the image . 

It may also be contemplated that on the occasion 
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of extracting the left end paint xl and the right end 
point x2 in the object area/ at the left end point 
extraction circuit 2202 acid at the right end point 
extraction circuit 2203,/ another condition for the 
secondary difference values SS(x) used in the 
extraction, for example such a condition that the 
secondary difference/ values SS(x) are not more than a 
fixed threshold, is added. 
( Embodiment 2 - 2 ) 

The present embodiment is arranged to apply, for 
example, the projection in the object area represented 
by Eq. (24) below to Eq. (21) in above Embodiment 2-1. 

... (24) 



f (x) = f c f(x,y)dy 

J b 



In this Eq. (24), "f(x,y) n represents the image 
data of the object area in the input image G, and "x" 
and "y" are coordinates on the horizontal axis and on 
the vertical axis, respectively. Further, "b" and "c" 
indicate the object area. 

When the projection in the object area represented 
by above Eq. (24) is used in above Eq. (21) in 
Embodiment 2-1, it can achieve the effects similar to 
those in the case where steps S2301 to S2304 described 
above are carried out in the averaged state of the data 
of the object area. Namely, mis judgment of presence or 
absence of the irradiation diaphragm can be prevented 
from being caused by alteration of coordinates of the 
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object area due "to influence of scattered rays, noise, 
and so on. Therefore, presence/absence of the 
irradiation diaphragm in the object area can be judged 
with better accuracy. 
5 (Embodiment 2-3) 

The present embodiment is applied, for example, to 
an image judgment device 2400 as illustrated in Fig. 
10. 

This image judgment device 2400 has the structure 

10 similar to that of the image judgment device 2100 in 

Embodiment 2-1 described above, but is different in the 
internal structure of the judgment unit 2110 therefrom. 

In the image judgment device 2400 of Fig. 10, 
elements operating similarly to those of the image 

15 judgment device 2100 of Fig. 7 are denoted by the same 
reference symbols and the detailed description thereof 
is omitted herein. Only the different structure from 
Embodiment 2-1 will be described in detail herein. 
The judgment unit 2140 is composed of, as 

20 illustrated in Fig. 10, a coordinate indication circuit 
2404, a secondary difference value calculation circuit 
2401 for calculating the sec/ondary difference values 
from data of an object area (either one of the areas A 
to D in the input image p of Fig. 8 in this example) in 

25 the input image according to coordinates indicated by 
the coordinate indication circuit 2404, a left end 
point extraction circuit 2402 for extracting a left end 

/ on 
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point of the irradiation area included in the object 
area, based on the secondary difference values 
calculated at the secondary difference value 
calculation circuit 2401, arid a right end point 
5 extraction circuit 2403 for extracting a right end 

point of the irradiation area included in the object 
area, based on the secondary difference values 
calculated at the secondary difference value 
calculation circui/c 2401. 
10 The judgment unit 2140 further has a memory 

circuit 2405 for storing coordinates of a plurality of 
left and right end points extracted for a plurality of 
lines by the left end point extraction circuit 2402 and 
the right end point extraction circuit 2403, an average 
15 coordinate calculation circuit 2406 for calculating an 
average coordinate at each end of the coordinates of 
the left and right end points stored in the memory 
circuit 2405, a first judgment circuit 2407 for judging 
whether the object area is an area with the irradiation 
20 diaphragm or an area without the irradiation area, 

using the averages of the coordinates of the left and 
right end points calculated at the average coordinate 
calculation circuit 2406, a coordinate variance 
calculation circuit 2408 for calculating variances of 
25 the coordinates of the left and right end points 

according to the judgment result of the first judgment 
circuit 2407, and a second judgment circuit 2409 for 
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again judging whether the object area is an area with 
the irradiation diaphragm or an area without the 
irradiation diaphragm, based on the variances 
calculated at the coordinate variance calculation 
circuit 2408. 

The coordinate indication Circuit 2404 is 
configured to give an indication of a coordinate for 
calculation of the secondary /difference value at the 
secondary difference value Calculation circuit 2401 
thereto after the left end point extraction circuit 
2402 and the right end opint extraction circuit 2403 
extract the left and r/ght end points. 

A processing program, for example, according to 
the flow chart as illustrated in Fig. 11, for 
controlling the operation of the judgment unit 2140 is 
preliminarily stored in the program memory 2130 and 
this processing program is read and executed by the 
control unit 2120, whereupon the judgment unit 2140 
operates as follows. 

First, supposing the area c/is an object area, the 
secondary difference value calculation circuit 2401 
calculates the secondary difference values SSi(x) from 
the data of the area C according to Eq. (25) below for 
each coordinate i indicates by the coordinate 
indication circuit 2404 /step S2501). 

SSi(x) = fi(x - d) - 2 x fi(x) + fi(x + d) 



(25) 
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In this Eq. (25), " fi(x)" represents image data of 
a one-dimensional line crossing the area C in the 
horizontal direction, "x" coordinates thereof, and 11 i" 
a coordinate of the horizontal line (a coordinate of 
the one-dimensional data) indicated by the coordinate 
indication circuit 2404. Further, "d" is a constant 
indicating a difference distance. 

Using the secondary difference values SSi(x) 
calculated at the secondary difference value 
calculation circuit 2403^, the left end point extraction 
circuit 2402 then ex-tracts the left end point xLi 
according to Eq. (^5) below (step S2502 ) . 

SSi(xLi) = min{SSi(x) 10 <, x <; x} 



x; coordinates on the horizontal axis between the 

left and right ends of the object area (coordinates 
between xa and xb) 

Then the memory circuit 2405 stores the left end 
point xLi extracted at the left end point extraction 
circuit 2402 (step S2503). 

Using the secondary difference values SSi(x) 
calculated at the secondary (difference value 
calculation circuit 2401, trie right end point 
extraction circuit 2403 aflso extracts the right end 
point xRi according to^Eq. (27) below (step S2504). 

SSi(xRi) = min{SSi(x) Ix <> x £ Length} 

... (27) 
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In this Eq. (27), "Length" indicates the length of 
the area C along the horizontal axis. 

Then the memory circuit 2405 stores the right end 
point xRi extracted at the right end point extraction 
circuit 2403 (step S2505 ) . 

After completion of the storage of the left end 
point xLi and the right end point xRi in the memory 
circuit 2405, the coordinate indication circuit 2404 
gives another instruction of a coordinate i of a new 
horizontal line to the secondary difference value 
calculation circuit 2401 (step S2506). 

According to this indication, the processing from 
step S2501 is carried out agai/. This loop processing 
is executed before the coordinate indication circuit 
2404 gives an indication o^ the end to the secondary 
difference value calculation circuit 2401. 

After steps S2501 to S2506 ar<s carried out 
repeatedly and the processing is/terminated according 
to the end indication from the /coordinate indication 
circuit 2404, the memory circuit 2405 is in a storage 
state of the left end points xLi and the right end 
points xRi corresponding /to the coordinates i indicated 
during the processing by the coordinate indication 
circuit 2404 to the Secondary difference value 
calculation circuit 2401 . 

Then the average coordinate calculation circuit 
2406 calculates the averages of the left end points xLi 
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and the right end points xRi stored in the memory 
circuit 2405 (step S2507). 

This obtains the average of the left end points 
xLi (left end point average) XL and the average of the 
5 right end points xRi (right end point average) XR. 

Next, similar to the irradiation diaphragm 
judgment circuit 2204 of Fig. 7, the first judgment 
circuit 2407 compares the left end point average XL and 
the right end point average XR calculated at the 

10 average coordinate calculation circuit 2406 with the 
left and right irradiation end points preliminarily 
obtained (the positions X a and X b of the irradiation 
ends on the horizontal axis H), and if the following 
relations are satisfied it judges that the left and 

15 right end points indicated by the left end point 

average XL and the right end point average XR are the 
left and right end points of the irradiation area, 
i.e., that the area C is an area with the irradiation 
area (or an area without the irradiation diaphragm). 

20 -e < XL - X a < e and 

-e < XR - X b < e 
On the other hand, if either one is not met, the 
judgment circuit judges that the area C is an area 
without the irradiation area (or an area with the 

25 irradiation diaphragm) (step S2508). 

The next steps S2509 and S2510 are carried out 
only if the first judgment circuit 2407 judges that the 

3^ 
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area is an area without: the irradiation diaphragm. 

Namely, the coordinate variance calculation 
circuit 2408 calculates the variances of the 
coordinates of the left and right end points, using the 
5 left end points xLi and the right end points xRi stored 
in the memory circuit 2405 and the left end point 
average XL and the right end point average XR 
calculated at the average coordinate calculation 
circuit 2406 (step S2509 ) . 
10 This operation provides the variance VL of the 

coordinates of the left end points and the variance VR 
of the coordinates of the right end points. 

The second judgment circuit 2409 compares the 
variances VL and VR of the coordinates of the left and 
15 right end points obtained at the coordinate variance 

calculation circuit 2408 with a predetermined threshold 
TH . When the following relations are satisfied, the 
second judgment circuit 2409 judges that the area C is 
an area with the irradiation area (or an area without 
20 the irradiation diaphragm). 
VL < TH and 
VR < TH 

On the other hand, if either one is not met, the second 
judgment circuit 2409 judges that the area C is an area 
25 without the irradiation area (or an area with the 
irradiation diaphragm) (step S2510). 

Steps S2501 to S2510 are also carried out for each 
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of the other areas A, B, D in the similar fashion to 
the area C. 

As described above, thes present embodiment is 
arranged to carry out the/detection of the irradiation 
ends on the plural linens in the object area while the 
coordinate indication circuit 2404 gives the 
indications of th^ coordinates i of one-dimensional 
data lines crossing the object area in the horizontal 
direction to/the secondary difference value calculation 
circuit 2^01. 

Since this structure is arranged to carry out the 
detection of the irradiation ends on the plural lines 
in the object area, the present embodiment is more 
unlikely to be affected by the scattered rays and the 
noise and can perform the judgment with better accuracy 
than in the case of the detection of the irradiation 
ends on a certain line. 

Since the variances of the irradiation ends are 
also used for the judgment of presence/absence of the 
irradiation diaphragm by provision of the second 
judgment circuit 2409, for example, when the image end 
portions are exposed to radiation because of the 
influence of the scattered rays, so as to result in 
making the irradiation ends extracted in the image end 
portions coincident with the irradiation ends 
preliminarily obtained, the degree of spread of the 
scattered rays can be quantified and thus mis judgment 
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of presence/absence of the irradiation diaphragm can be 
prevented for sure. 
( Embodiment 2 - 4 ) 

The present embodiment is Arranged to calculate 
5 the secondary difference values SS(x) according to 

above Eq. (21) and to extract the left end point xl of 
the irradiation area included in the object area, using 
the secondary differences values SS(x), for example, 
similar to Embodiment /2-1 described above. At this 
10 time the present embodiment also uses the sign of the 

primary difference/Value S(xl) given by Eq. (28) below. 

S(xl) = f(xl) - f(xl - d) ... (28) 

Specifically, for example, if the sign of the 
primary difference value S(xl) is "negative" and if 
15 above Eq. (22) is satisfied, the left end point xl is 
regarded as a left end point of the irradiation area. 

For extracting the right end point x2 of the 
irradiation area included in the object area, the 
present embodiment also uses the sign of the primary 
20 difference value S(x2) given by Eq. (29) below. 

S(x2) = f(x2 + d) - f(x2) ... (29) 

For example, if the sign of the primary difference 
value S(x2) is "negative" and if above Eq. (23) is 
satisfied, the right end point x2 is regarded as a 
25 right end point of the irradiation area. 

Using the left end point xl and the right end 
point x2 obtained in this way, the judgment of 

35 
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presence/absence of "the irradiation diaphragm in the 
object area is carried out in the similar fashion to 
that in Embodiment 2-1 described above. 

When the conditions of the signs of the primary 
5 difference values are also added to the detection of 
the irradiation ends of the object area, 
p reS ence/absence of the irradiation diaphragm in the 
object area can be judged with consideration to the 
inclination of image data outside the irradiation area 
10 due to the scattered rays. This permits the 

presence/absence of the irradiation diaphragm in the 
object area to be judged without mis judgment and with 
better accuracy even if there is a portion or the like 
where densities change quickly in the object area. 
15 ( Embodiment 2 - 5 ) 

The present embodiment i^s arranged first to 
calculate the secondary difference values from the data 
of the object area, fpr example similar to Embodiment 
2-1, but at this tilne, the present embodiment is 
20 arranged to calculate the secondary difference values 

from data o>tained after the data of the object area is 
f ilterec 

Specifically, for example, the image data of a 
one-dimensional line of an object area is represented 
25 by "f(x) n , the data is subjected to a filtering process 



according to Eqs. (30) and ('31) below, and the 
secondary difference value's are calculated from the 
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values F2 obtained/ as a result. 



Fl(x) = min{f(x + xl ) - h(xl)|-d < xl < d} 



(30) 



F2(x) = max{Fl(x - xl ) + h(xl)|-d < xl < d} 




(31) 



-d z x ^ d 
otherwise 



15 



20 



When Embodiment 2-1 is modified so as to calculate 
the second difference values after the data of the 
object area is smoothed by the filtering process as 
described above, the presence/absence of the 
irradiation diaphragm in the object area can be judged 
with better accuracy without influence of the noise, 
particularly, without influence of the noise on the 
line. 

( Embodiment 3 - 1 ) 

The present embodiment is carried out, for 
example, by an image judgment device 3100 as 
illustrated in Fig. 12. 

This image judgment device 3100 is composed of a 
judgment unit 3110, a control unit 3120 for controlling 
the operation of the judgment unit 3110, and a program 
memory 3130 to which the control unit 3120 makes 
access, as illustrated in Fig. 12. 

The judgment unit 3110 is composed of a coordinate 
indication circuit 3201, a characteristic value 
extraction circuit 3202 for extracting a characteristic - 
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value from data of an object area in an input image 
according to a coordinate indicated by the coordinate 
indication circuit 3201, an end point extraction 
circuit 3203 for extracting a coordinate of an end 
5 point of the irradiation area included in the object 
area, based on the characteristic value extracted at 
the characteristic value extraction circuit 3202, a 
memory circuit 3204 for storing coordinates of end 
points extracted at the end point extraction circuit 

10 3203, and a judgment circuit 3205 for judging whether 
the object area is an area with the irradiation 
diaphragm or an area without the irradiation diaphragm, 
from the coordinates of the end points stored in the 
memory circuit 3204. 

15 Here, for example, an image G illustrated in Fig. 

13 is supplied as an input image to the judgment unit 
3110. This input image G is a two-dimensional, 
radiographic, thoracic part image obtained by 
photographing the thoracic part with the irradiation 

20 diaphragm by use of an imaging device with the 
irradiation diaphragm function. 

In above Fig. 13, "X a " and "X b " represent positions 
of irradiation area ends on the horizontal axis X. "A" 
to "D M are areas at the edges of the image part. In 

25 this example, the areas A and B out of these "A" to "D" 
are areas with the irradiation diaphragm while the 
other areas C and D are areas without the irradiation 
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diaphragm. 

Various processing programs for controlling the 
operation of the judgment unit 3110 are preliminarily 
stored in the program memory 3130. 

Specifically, for example, a processing program 
according to the flow chart as illustrated in Fig. 14 
is preliminarily stored in the program memory 3130 and 
this processing program is read and executed by the 
control unit 3120 to operate the judgment unit 3110 as 
follows . 

Described below is a process of judging whether 
the irradiation diaphragm is present or absent at the 
lower edge . 

First, the coordinate indication circuit 3201 
gives indications of plural coordinates on the X-axis 
between X a and X b in the input image G of above Fig. 13 
to the characteristic value extraction circuit 3202 
(step S3301). 

Specifically, the coordinate indication circuit 
3201 successively indicates, for example, coordinates i 
(= 1 to 10) on the X-axis, which are ten points 
distributed at equal intervals between X a and X b . 

According to each coordinates i indicated by the 
coordinate indication circuit 3*201, the characteristic 
value extraction circuit 3202/then calculates a 
characteristic value of one/dimensional image data 
corresponding to the coordinate i, for example, the 
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secondary difference values SSi(y) according to Eq. 
(41) below ( step /S3302 ) . 

SSi(y) = fi(y - d) - 2 x fi(y) + fi(y + d) 

... (41) 

In this Eq, (41), " fi(y)" indicates image data of 
a one-dimensional line crossing the object area in the 
vertical direction, "y" coordinates thereof, and "d" a 
constant indicating a difference distance. 

Using the secondary difference values SSi(y) 
obtained at the characteristic value extraction circuit 
3202, the end point/extraction circuit 3203 then 
extracts a coordinate yi of an end point ( an 
irradiatioi/end point ) of the irradiation area included 
in the x^bject area according to Eq. (42) below (step 
S3303 ) . 

SS(yi) = min{SSi(y)|0 < y < Length} 

... (42) 

In this Eq. (42), "Length" represents the length 
of the input image G along the vertical axis ( the 
direction of the Y-axis). 

Then the memory circuit 3204 stores the coordinate 
yi of the irradiation end point obtained at the end 
point extraction circuit 3203 (step S3304). 

The above processing of steps S3301 to S3304 is 
carried out repeatedly before the coordinate i 
indicated by the coordinate indication circuit 3201 
reaches "10". This results in storing the coordinates 
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yl, y2, y3,..., ylO of the irradiation end points 
corresponding to the coordinates i = 1, 2, 3,..., 10 
indicated by the coordinate indication circuit 3201, in 
the memory circuit 3204. 
5 After completion of the processing of steps S3301 

to S3304 up to the coordinate i = 10, the judgment 
circuit 3205 then calculates an average of the 
coordinates yl, y2, y3,...> ylO stored in the memory 
circuit 3204 (step S3305) and calculates a variance Bv 

10 of the irradiation end points (step S3306 ) . 

Then the judgment circuit 3205 judges that the 
object area is an area without the irradiation 
diaphragm, if the variance Bv is not less than a 
predetermined threshold TH; otherwise, it is judged 

15 that the object area is an area with the irradiation 
diaphragm (step S3307 ) . 

The threshold TH is a constant which is determined 
experimentally. 

As described above, the present embodiment is 

20 arranged to carry out the extraction of irradiation 
ends of plural column lines in the object area for 
which presence/absence of the irradiation diaphragm is 
judged and judge the presence/ absence of the 
irradiation diaphragm of the object area according to 

25 the variance of the irradiation ends. If there is an 
irradiation area in the object area the coordinates of 
the irradiation end points of the plural column lines 
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will be aligned approximately on one horizontal axis 
and the variance thereof will be thus small. If there 
is no irradiation area in the object area the 
coordinates of the irradiation end points of the plural 
5 column lines will be distributed and the variance will 
be thus large. 

Therefore, the presence/absence of the irradiation 
diaphragm in the object area can be determined with 
accuracy because of the structure using the variance. 
10 Since the device is constructed to iase the 

secondary difference value in order to /extract the 
irradiation end, boundary points can JzSe extracted with 
accuracy between an area irradiated yairectly and the 
other areas even in a photographic /image obtained by 
15 photographing a subject with low ycransmittances of 
radiation. Therefore, the device of the present 
embodiment can judge the presemce/absence of the 
irradiation diaphragm in the object area including the 
irradiation area with accuracy. In addition, the 
20 presence/ absence of the irradiation diaphragm in the 
object area can be judged /with accuracy even in a 
photographic image in whyLch a portion with low 
radiation transmittance's such as the abdominal part or 
the like overlaps witn an end portion of the image. 
25 Although the present embodiment is arranged to use 

the secondary difference value i^n order to make a 
judgment of the irradiation diaphragm, the apparatus of 
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the present: invention does noy always have to be 
limited to this; for example/ where change of density 
is quick at an irradiation/end, the apparatus may also 
be arranged to use the pi/imary difference values or 
higher-order differences values. In this case, the 
primary difference values or the higher-order 
difference values are obtained from the object area and 
a first appearing/ point of a value not less than a 
predetermined threshold is regarded as a candidate for 
an irradiation end ♦ 

Although the present embodiment is arranged to use 
the variance in order to make a judgment of the 
irradiation diaphragm, the apparatus of the present 
invention does not always have to be limited to this; 
for example, the apparatus may also be arranged to use 
another index that can indicate discrete degrees of 
frequency, such as standard deviation. 
(Embodiment 3-2) 

The present embodiment is arranged to apply the 
projection of a constant width in an object area 
according to Eq. (43) below to Eq. (41) in Embodiment 
3-1 described above. 

... (43) 



fi(y) = / f(x,y)dx 

Jxi-l 



In this Eq. (43), "fi(y)" represents image data of 
a one-dimensional line in the object area in above Eq. 
(41), and "Xi" and "Xi-1" represent coordinates 
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indicated by the coordinate indication circuit 3201. 

When Embodiment 3-1 is modified so as to apply the 
projection in the object area represented by above Eq. 
(43) to above Eq. (41) as described above, the 
5 modification can enjoy the same effects as in the case 
where aforementioned steps S3301 to S3307 are carried 
out in an averaged state of the data of the object 
area. Namely, this structure can prevent mis judgment 
about the presence/absence of the irradiation diaphragm 
10 from being caused by change of coordinates of the 

object area due to influence of the scattered rays, the 
noise, and so on. The presence/absence of the 
irradiation diaphragm in the object area can be judged 
with better accuracy accordingly. 
15 (Embodiment 3-3) 

The present embodiment is arranged to calculate 
the secondary difference value/ SSi(y) according to 
above Eq. (41) and extract the irradiation end of the 
object area using the secondary difference values 
SSi(y), for example, as /Embodiment 3-1 described above 
was. At this time, ttte present embodiment also uses 
the sign of the primary difference value Si(y) 
expressed by Eq. y( 44 ) below. 

Si(y) = fi(y) - fi(y - d) ... (44) 

Specifically, for example, if the sign of the 
primary difference value Si(y) is "negative" and if 
above Eq. (42) is satisfied, the point "y" is regarded 
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as an irradiation end. 

Using the irradiation ends obtained in this way, 
whether the irradiation diaphragm is present or absent 
in the object area is judged in the similar fashion to 
aforementioned Embodiment 3-1. 

When Embodiment 3-1 is modified to use the sign of 
the primary difference value as well on the occasion of 
detecting the irradiation end in the object area, the 
presence/absence of the irradiation diaphragm in the 
object area can be judged, also taking the inclination 
of image data outside the irradiation area due to the 
scattered rays into consideration. Therefore, whether 
the irradiation diaphragm is present or absent in the 
object area can be judged with better accuracy without 
misjudgment, even if there is a quickly changing 
portion of density or the like within the object area. 
( Embodiment 3 - 4 ) 

The present embodiment is ^arranged first to 
calculate the secondary difference values from the data 
of the object area, similar to Embodiment 3-1; but at 
this time, the data o£^the object area is subjected to 
a filtering process and the secondary difference values 
are calculated/from the data after the filtering 
process, 

Specifically, for example, where the image data of 
a one-dimensional line of the object area is "f(x) H , it 
is subjected to the filtering p/ocess according to Eq. 
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(45) and Eq. (46) b^low and the secondary difference 
values are calculated from values F2 obtained as a 
result. / 

Fli(y) = min{fi(y + yl ) - h(yl)|-d < yl < d} 

... (45) 

F2i(y) = max{Fli(y - yl ) + h(yl)|-d < yl < d} 

... (46) 

h (v\ = r 0 : -d ^ x ^ d 
nix; 1 -oo : otherwise 

When Embodiment 3-1 is modafied so as to calculate 
the secondary difference values after the data of the 
object area is smootheci^y the filtering process as 
described above, wh^her the irradiation diaphragm is 
present or absept in the object area can be judged with 
better accuracy without being affected by the noise, 
particularly , without being affected by the noise on 
the l;me. 
( Embodiment 4-1 ) 

Fig. 15 is a block diagram to show the structure 
of an angle extraction device and an area extraction 
device of two-dimensional image data according to 
Embodiment 4-1. In Fig. 15, reference numeral 4101 
designates a characteristic quantity calculation means 
for calculating a characteristic quantity from data of 
a one-dimensional line at a coordinate indicated by a 
coordinate indication means 4103, 4102 an end point 
extraction means for extracting an area end portion 
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from the characteristic quantity calculated at the 
characteristic quantity calculation means 4101, 4103 a 
coordinate indication means for indicating a coordinate 
of image data to be calculated by the characteristic 
5 quantity calculation means 4101, 4104 an end point 

memory means for storing a coordinate of an end point 
extracted at the end point extraction means 4102, and 
4105 a rotation angle indication means for indicating a 
rotation angle of a rotation axis onto which 

10 coordinates of end points stored in the end point 
memory means 4104 are projected. 

Reference numeral 4106 denotes an accumulated 
quantity calculation means for calculating an 
accumulated quantity of the coordinates of the end 

15 points projected onto the aforementioned rotation axis 
and stored in the end point memory means 4104, 4107 an 
accumulated quantity memory means for storing the 
accumulated quantity calculated at the accumulated 
quantity calculation means 4106, 4108 a rotation angle 

20 judgment means for judging a rotation angle from the 

accumulated quantity stored in the accumulated quantity 
memory means 4107, and 4109 an end point judgment means 
for extracting an area end from the angle judged at the 
rotation angle judgment means 4108 and the accumulated 

25 quantity stored in the accumulated quantity memory 
means 4107. 

Fig. 16 is a flow chart of a processing sequence 
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in the angle extraction device and the area extraction 
device in Embodiment 4-1. Fig. 17 is a schematic 
diagram for explaining the processing in the angle 
extraction device and area extraction device of the 
present embodiment, in which "a" represents two- 
dimensional data (data within a light receiving surface 
in the case of copiers etc., or data in a sensor 
surface in the case of X-ray apparatus etc.), "b" a 
rectangular area (an irradiation area in the case of 
the X-ray apparatus or a sheet surface in the case of 
the copiers etc.), "c H a start point of the rotation 
axis, "d" an area end point, "e" a start axis of angle, 
" f" a rotation axis, and "g" an angle difference Gn 
between the start axis e of angle and the rotation axis 
f . 

The operation of the present embodiment will be 
described below according to the flow chart of Fig. 16. 

The characteristic quantity calculation means 4101 
of Fig. 15 calculates the secondary difference values 
according to Eq. (51) below (steps S4101 and S4102). 
In the equation fi(x) represents image data of the i-th 
row line indicated by the coordinate indication means 
4103, SSi(x) the secondary difference values, and "c" a 
npastant . 

SSi(x) = fi(x - c) - 2 x fi(x) + fi(x + c) 

... (51) 

Then the end point extraction means 4102 extracts 

41 
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a candidate for an end point of the area according to 
Eq. (52) below (step S4103 ) . 

For the left area end of the i-th row, the 
following equation is used. 
5 xi = min{SSi (x) I 0 <; x <; x} • * * ( 52 ) 

In this equation x indicates a coordinate on an axis 
along an indicated direction within the area (the 
details will be described hereinafter and reference 
should be made to Eq. (76)). Then the end point memory 

10 means 4104 (Fig. 15) stores the coordinate of the area 
end obtained according to Eq. (52) (step S4104). Then 
steps S4101 to S4104 are repeated before the processing 
is complete for all rows indicated by the coordinate 
indication means 4103 (Fig. 15). 

15 Next, the rotation angle indication means 4105 

(Fig. 15) indicates the rotation angle Gn (Fig. 17) 
(step S4105), and the accumulated quantity calculation 
means 4106 (Fig. 15) performs processing indicated by 
Eq. (53) to Eq. (57) below to calculate the accumulated 

20 quantity L(Gn,X) of coordinates x on the aforementioned 
rotation axis f (Fig. 17) (step S4106). The projection 
of the area end points onto the aforementioned rotation 
axis is obtained from projected values of candidate 
points onto the rotation axis rotated about the start 

25 point c(xc,yc) of the aforementioned rotation axis. 

x ±en = "( xi - xc)*cos(6n) + (Yi - yc)*sin(6n) 

... (53) 
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Here, Xi, Yi represent coordinates of a candidate point 
for the area end of the i-th row, and X i6n represents 
projected positions of the candidate points for the 
area ends on the rotation axis x, which are points 
resulting from rotation by 6n about the start point c 
(Fig. 17) of the aforementioned rotation axis. 

m 

L(0n,x) =£)Z(X / X i8n ) (54) 
i=l 

Here, L(0n,X) indicates the number of projected 
candidate points on the axis resulting from the 
rotation of the aforementioned rotation axis by 6n. 

Here, Z(X,Xi) =1, X - dl < Xi < X + dl 

= 0, otherwise, ... (55) 

where dl is a constant. 

Then the accumulated quantity memory means 4107 
(Fig. 15) stores the accumulated quantity L(0n,x) 
obtained by Eq. (54) (step S4107 ) . Further, the 
processing of steps S4105 to S4107 is repeated for all 
angles 6n indicated by the rotation angle indication 
means 4105 (Fig. 15). 

Then the rotation angle judgment means carries out 
the processing of Eq. (56) to Eq. (58) below to extract 
the rotation angle of the area (step S4108). 

n(9n) = L(Gn,Xmax - 1) + L(0n,Xmax) + L(9n,Xmax + 1) 

... (56) 

This represents the number of overlaps of area ends at 
the maximum overlap point in the projection of area 
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ends onto the rotation axis. Here, Xmax indicates a 
coordinate of X that satisfies the following, 

L(Gn,Xmax) = max{L( 9n, X ) | X6kl} ... (57) 

Here, kl represents coordinates on the rotation axis. 
5 After all, the rotation angle 6 is an angle 

obtained according to the following equation. 

9 = max{n(6n) | GnGK} ... (58) 

In this equation K is an arbitrary domain of 
definition. 

10 Then the end point judgment means 4109 (Fig. 15) 

extracts Xmax satisfying Eq. (59) as a candidate for a 
left end point of the rectangular area from 0 
determined by Eq. (58) and the accumulated quantities 
L(6,X) stored in the accumulated quantity memory means 

15 4107 (Step S4109). 

L(9,Xmax) = max{ L( 6 , X ) | XGkl } ... (59) 
Likewise, candidates are also extracted for the 
right, upper, and lower end points. It is, however, 
noted that they do not always have to be extracted, 

20 because the angle of the area is already determined . 

As described above, the present embodiment has the 
following effect because of the use of the secondary 
difference value; when the area is extracted, boundary 
points between the area directly irradiated with X-rays 

25 and the other areas can be extracted as end points of 
the area with accuracy even from an object with low X- 
ray transmittances . 
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Further, since the angle and end points of the 
area are determined from the accumulation of end 
points, the angle and end points of the area can be 
extracted with accuracy and within short computation 
5 time. 

In the case of the copiers, FAX, OCR, etc., a 
necessary area can be extracted based on the area and 
angle obtained. Further, the present embodiment has 
the effect of capability of automatically carrying out 
10 selection of a sheet with accuracy, because the size of 
the area can be identified. In addition, even if a 
sheet is placed obliquely, correction can be made and 
thus the result is equivalent to that where the sheet 
is placed at a normal position. 
15 (Embodiment 4-2) 

In the present embodiment; where the image data of 
a one-dimensional line is/£(x) and the secondary 
difference values the??^of are SS(x) defined by Eq. (1), 
the sign of the primary difference value S(XL) of Eq. 
20 (60) below is/also added to the operation of extracting 
the end opint XL at the end point extraction means 4102 
(Fig^Xs). 

S(XL) = f(XL) - f(XL - d) ... (60) 

For example, for extracting the left end point, XL 
25 that makes S(XL) negative and that satisfies Eq. (61) 
below is regarded as a left end point. 

SS(XL) = min{SS (x) I 0 <; x <; x} (61) 
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As described above, the present embodiment can 
also take the inclination of the image data outside the 
area into consideration, because the sign of the 
primary difference is added to the extraction of the 
area end. This eliminates erroneous extraction of a 
quickly changing portion of density in the area and 
thus the area end point can be extracted with better 
accuracy . 
( Embodiment 4-3 ) 

In the present embodiment, whefe the one- 
dimensional data of the i-th row i:s fi(x) and is 
subjected to the filtering process according to Eqs . 
(64) and (65) below and valuesr resulting from the 
filtering process are defined as Fli(x) and F2i(x), the 
characteristic quantity calculation means 4101 (Fig. 
15) uses F2i(x) for calculating the secondary 
difference values defi/ned by Eq. (51). 

Fli(x) = min{fi(x + xl ) - h(xl)|-d < xl < d} 

... (64) 

Further, F2i(x) is defined as follows. 

F2i(x) = max{Fli(x - xl ) + h(xl)|-d < xl < d} 

... (65) 

Here, h(x) is a function defined below. 
h(x) =0, -d < xl < d 

= oo, otherwise ... (66) 

Hence, the secondary difference values are 
calculated as folio 
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SSi(x) = F2i(x - c) - 2 x F2i(x) + F2i(x + c) 

... (67) 

As described above, the/present embodiment is 
arranged to smooth the on^-dimensional image data for 
calculation of the secondary difference values by the 
filtering process, ylfhereby accomplishing the effect of 
being not affected by the noise, particularly by the 
noise on the y&ine. 
(Embodiment 4-4) 

The present embodiment is an example in which the 
start point c (Fig. 17) of the aforementioned rotation 
axis is placed at the barycenter of data not less than 
a fixed density. The flow of processing will be 
described according to a flow chart of a procedure 
sequence for obtaining an approximate center in the 
area, shown in Fig. 18. 

First, an accumulation histogram of whole image 
data is produced in step S4201 and upper densities 
above TH % are extracted (steps S4202, S4203, S4204 ) . 
Then the barycenter of the upper densities above TH % 
is regarded as an approximate barycenter of the 
rectangular area. 

- , f ff (x,y)xdxdy ... (68) 

JJf <x,y)dxdy 

- = JJf (x /y )ydxdy ... (69) 
y JJf <x,y) dxdy 
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Here, f(x,y) represents an image of the upper 
densities above TH % and y, x stand for the barycenter. 

As described above, since the present embodiment 
is arranged to calculate the barycenter of data above 
5 the fixed density value, the present embodiment can 
extract the approximate barycenter of the area of 
interest with accuracy without being affected by the 
peripheral data, even if the object area is in the 
peripheral part of the two-dimensional image data. 

10 Further, since the start point c (Fig. 17) of the 

aforementioned rotation axis f is placed in the object 
area, degrees of change in accumulated quantities of 
area ends become large against rotation of the rotation 
axis and the present embodiment can thus enjoy the 

15 effect of capability of extracting the rotation angle 
and area with better accuracy. 
(Embodiment 4-5) 

The present embodiment is arranged as follows. 
When the end point extraction means 4102 (Fig. 15) 

20 extracts an end point, a point is not extracted as a 

candidate for an end point if f maxi < Th at the point, 
f maxi being defined below. 

f maxi = max{f i(x) | x€k} ... (70) 
Here, "k" indicates a domain of definition of row data 

25 and Th a fixed threshold value. 

As described above, the present embodiment can 
eliminate end points that can be points outside the 
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object: area with high possibilities and thus can 
extract the area and the angle of the area with better 
accuracy. 
(Embodiment 4-6) 
5 The present embodiment uses a first order 

differential or a higher order differential as a 
characteristic quantity at the characteristic quantity 
calculation means 4101 (Fig. 15). When Di(x) 
represents differential values of the i-th row, the end 
10 point extraction means 4102 (Fig. 15) extracts a point 
satisfying the condition of xi > th as a candidate for 
an end point on the occasion of extracting a candidate 
for the left end point of the i-th row. 

As described above, the present embodiment has the 
15 effect of shorter computation time, because it adopts 
the simple differentiation. 
(Embodiment 4-7) 

Fig. 19 is a block diagram to /show the structure 
of an area extraction device according to the present 
20 embodiment. In Fig. 19, reference numeral 4301 

designates a secondary difference value calculation 
means for calculating the secondary difference values 
of one-dimensional image data in a designated direction 
(for example, which is determined according to the 
25 value determined by th^ angle extraction device 

described above), 43tf2 a left end point extraction 
means for extracting a left end point of the area, 
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based on the secondary difference values calculated at 
the secondary difference/value calculation means 4301, 
~7 and 4303 a right end rioint extraction means for 

C ^/ extracting a right y^nd point of the area, based on the 

5 secondary difference values calculated at the secondary 
difference value calculation means 4301. Fig. 20 is a 
flow chart or a processing procedure sequence in the 
area extraction device according to the present 
embodinjent . 

10 The flow of the processing in the present 

embodiment will be described according to Fig. 20. The 
secondary difference value calculat/on means 4301 (Fig. 
^ 19) calculates the secondary difference values SS(x) 
C^l / according to a calculation equation defined by Eq. (71) 

15 below (step S4301 ) . Here, fjfa) represents the one- 
dimensional data of a line/crossing the area in the 
designated direction ancr x represents coordinates 
thereof. Further, "<y denotes a constant indicating a 
difference distance 
20 SS(d) = f(x - d) - 2 x f(x) + f(x + d) 

... (71) 

The left end point extraction means 4302 (Fig. 19) 
extracts the left end point xl according to Eq. (72) 
below (step S4302). Here, x represents a coordinate on 

25 the horizontal axis in the area. 

SS (xl) = min{SS(x) I 0 <. x <; x} • * - (?2) 

Then the right end point extraction means 4203 
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(Fig. 19) extracts the right end point x2 according to 
Eq. (73) below (step S4303 ) . 

SS(x2) = min{SS(x) Ix <s x £ Length} 



Here, "Length" indicates the length of the image data 
along the horizontal axis. 

As described above, s/nce the present embodiment 
uses the secondary difference values, it has the effect 
of capability of extracting the area with accuracy even 
if the density changer is gentle at the area end. 
(Embodiment 4-8) 

The present embodiment uses the projection of the 
image area as defined by Eq. (74) below, as the one- 
dimensional data f(x) in aforementioned Eq. (71). 



Here, f(x,y) represents the image data in the 
designated direction and x, y coordinates on the 
horizontal and vertical axes, respectively (which are 
inclined in a certain direction). Further, "b", "c" 
represent an arbitrary domain across the object area. 

As described above, since the present embodiment 
uses the projection in the image area, it has the same 
effect as in the case of averaging the data of the 
image area. Therefore, the present embodiment can 
prevent the extraction from being affected by the noise 
inside and outside the object area and thus can extract 



(73) 




(74) 
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■the area with better accuracy. 
( Embodiment 4-9 ) 

Fig. 21 is a block diagram to show the /structure 
of an area extraction device according to /the present 
5 embodiment. In Fig. 21, reference numeral 4401 

designates a secondary difference value /calculation 
means for calculating the secondary difference values 
of one-dimensional image data in a designated direction 
(for example, which is determined according to the 

10 value determined by the angle extraction device 

described above) indicated by a coordinate indication 
means 4404, 4402 a left end point/extraction means for 
extracting a left end point of tfre area, based on the 
secondary difference values calculated at the secondary 

15 difference value calculation mefans 4401, a right end 

point extraction means for ext/racting a right end point 
of the area, based on the seaondary difference values 
calculated at the secondary /difference value 
calculation means 4401, and 4404 the coordinate 

20 indication means for indicating a coordinate of one- 
dimensional data for calculation of the secondary 
difference values at the secondary difference value 
calculation means 4401 /after extraction of the right 
area end at the right/end point extraction means 4403. 

25 Numeral 4405 denotes a memory means for storing 

coordinates of right end points and left end points 
extracted at the left end point extraction means 4402 
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and at -the right; end point extraction means 4403, and 
4406 an average coordinate calculation means for 
calculating an average coordinate for each end point 
from the left and right ends stored in the memory means 
4405. 

Fig. 22 is a flow chart of a processing procedure 
sequence of the area extraction device according to the 
present embodiment. The operation will be described 
according to Fig. 22. 

The secondary difference value ^calculation means 
4401 (Fig. 21) calculates the secondary difference 
values SSi(x) according to a calculation equation 
defined by Eq. (75) below ( stej/ S4401 ) . Here, fi(x) 
represents one-dimensional darca of a line crossing the 
area in the designated direction, x coordinates 
thereof, and 11 i" a coordinate of the row in the 
designated direction indicated by the coordinate 
indication means 4404/ Here, "d" represents a constant 
indicating a difference distance. 

SSi(x) = fi(x - d) - 2 x fi(x) + fi(x + d) 



Then the left end point extraction means 4402 
(Fig. 21) extracts the left area end XLi according to 
Eq. (76) below (step S4402 ) . 

Here, x represents a coordinate on the axis along 

the designated direction in the area. 



(75) 
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SSi (XLi) = min{SSi (x) I 0 <; x ^ x} 

... (76) 

Then the memory means 4405 (Fig. 21) stores the 
left area end XLi extracted at the left end point 
5 extraction means 4402 (step S4403 ) . 

Next, the right end point extraction means (Fig. 
21) extracts the right area end XRi according to Eq. 
(77) below (step S4404 ) . 

SSi (XRi) = min{SSi (x) Ix <> x ^ Length} 
10 ... (77) 

Here, "Length" represents the length of the image 
data along the horizontal axis in the predetermined 
direction. 

Then the memory means 4405 (Fig. 21) stores the 
15 right area end XRi extracted at the right end point 
extraction means 4403 (step S4404). 

Next, the coordinate indication means 4404 (Fig. 
21 ) indicates a coordinate of a row in a new designated 
direction and step S4401 to step S4405 are repeated. 
20 An indication of the end of the loop is issued from the 
coordinate indication means 4404 (Fig. 21) (step 
S4406) . 

Next, after completion of the above extraction of 
the left and right end points, the average coordinate 
25 calculation means 4406 (Fig. 21) calculates an average 
of the coordinates of the left and right end points 
stored in the memory means 4405 ( step S4407 ) . 
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As described above, jsince "the present embodiment 
uses the secondary difference values, the boundary 
points of the area /Can be extracted with accuracy even 
if the density values vary gently at the boundary of 
the area. f 

Further, since the area end is determined from the 
average of the area ends on plural lines, the present 
embodiment has the effects of being more resistant to 
the noise and capable of extracting the area with 
higher accuracy than in the case using a single point. 
(Embodiment 4-10) 

In the present embodiment, wbiere the image data of 
a one-dimensional line in the designated direction is 
expressed by f(x) and the secondary difference values 
thereof by SS(x) defined by/Eq. (78) below, the sign of 
the primary difference va^ue S(XL) defined by Eq. (79) 
below is added to the extraction of the left end point 
XL at the left end po/nt extraction means 4302 (Fig. 
19) or 4402 (Fig. 2JL\ ) . 

SS(x) = f(x - d) - 2 x f(x) + f(x + d) 

... (78) 

S(XL) = f(XL) - f(XL - d) ... (79) 

For example, for extracting a left end point, XL 

that makes s(XL) negative and that satisfies Eq. (80) 

below is regarded as a left end point. 

SS(XL) = min{SS (x) I 0 <; x <; x} (80) 
S(XR) = f(XR + d) - f(XR) ... (81) 
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Likewise, where the right end point is extracted 
at the right end point extraction means 4303 (Fig. 19) 
or 4403 (Fig. 21), XR that makes S(XR) defined by above 
Eq. (81) negative and that satisfies Eq. (82) below is 
5 regarded as a right area end. 

SS (XR) = min{SS(x) Ix ^ x £ Length} 

... (82) 

As described above, since the present embodiment 
is arranged to add the sign of the primary difference 
10 to the extraction of the area end, it can also take the 
inclination of the image data outside the area into 
consideration. Therefore, erroneous extraction can be 
prevented at a quickly changing portion of density in 
the area, so that the area can be extracted with better 
15 accuracy. 

(Embodiment 4-11) 

In the present embodiment , /where the one- 
dimensional data in the designated direction is defined 
by f(x) and values resulting from the filtering process 
20 according to Eqs. (83), /84) below are defined by 

Fl(x), F2(x), the values F2(x) are used for calculation 
of the secondary difference values at the secondary 
difference value/calculation means 4301 (Fig. 19) or 
4401 (Fig. 21 ), 
25 Fl(x) = min{f(x + xl ) - h(xl)|-d < xl < d} 

... (83) 

Further, F2(x) is defined as follows. 
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F2(x) = max{Fl(x - xl ) + h( xl ) | -d < xl < d} 

... (84) 

Here, h(x) is a function defined as follows. 
h(x) = 0, -d < xl < d 

= -«>, otherwise ... (85) 

As described above ,/since the present embodiment 
is arranged to smooth/the one-dimensional image data 
for the computation of the secondary difference values 
by the filtering process, it has the effect of being 
10 not affected by the noise, particularly, by the noise 
on the Lme . 

While the present invention has been described 
with respect to what is presently considered to be the 
preferred embodiments, it is to be understood that the 
15 invention is not limited to the disclosed embodiments. 
To the contrary, the invention is intended to cover 
various modifications and equivalent arrangements 
included within the spirit and scope of the appended 
claims. The scope of the following claims is to be 
20 accorded the broadest interpretation so as to encompass 
all such modifications and equivalent structures and 
functions . 



